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Abstract 

In this paper, we study a boundary control problem associated to the stationary Rayleigh-Benard- 
Marangoni (RBM) system in presence of controls for the velocity and the temperature on parts of the 
boundary. We analyze the existence, uniqueness and regularity of weak solutions for the stationary 
RBM system in polyhedral domains of R®, and then, we prove the existence of the optimal solution. 
By using the Theorem of Lagrange multipliers, we derive an optimality system. We also give a second- 
order sufficient optimality condition and we establish a result of uniqueness of the optimal solution. 
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1 Introduction 

Fluid movement by temperature gradients, also called thermal convection, is an important process in 
nature. Its main applications appear in industry, as for instance, in the growth of semiconductor crystals, 
but also, thermal convection is the basis for the interpretation of several natural phenomena such as the 
movement of the earth’s plates, the solar activity, large scale circulations of the oceans, movement in the 
atmosphere, among others. A model of particular interest consists of a horizontal layer of a fluid in a 
container heated uniformly from below, with the bottom surface and the lateral walls rigid and the upper 
surface open to the atmosphere. Due to heating, the fluid in the bottom surface expands and it becomes 
lighter than the fluid in the upper surface, so that, by effect of the buoyancy, the liquid is potentially 
unstable. Because of the instability, the fluid tends to redistribute. However, this natural tendency will 
be controlled by its own viscosity. On the other hand, the upper surface, which is free to the atmosphere, 
experiences changes in its surface tension as a result of the temperature gradients in the surface. Then, it 
is expected that the temperature gradient exceeds a critical value, above which the instability can manifest. 

The first experiments to demonstrate the beginning of Thermal instability in fluids were developed by 
Henri Benard in 1900 (see [5]). In his experiments, Benard considered a very thin layer of liquid, around 
1 mm of depth, in a metal plate maintained at a constant temperature. The upper surface was usually 
free and it was in contact with the air, which was at a lower temperature. Benard experimented with a 
variety of liquids with different physical characteristics, mainly interested in the effect of viscosity on the 
convection, using liquids of high viscosity like wax whale melted and paraffin. In all these cases, Benard 
found that when the temperature of the plate gradually increased, at a certain moment, the layer lost 
stability and formed patterns of hexagonal cells, all alike and correctly aligned. 
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A first theoretical interpretation of thermal convection was provided by Lord Rayleigh in 1916 (see 
[37]), whose analysis was inspired by Benard’s experiment. Rayleigh assumed that the fluid was confined 
between two horizontal thermally conductive plates and the fluid was being heated from below. Rayleigh 
considered that the effect of buoyancy is the only one responsible of the beginning of the instability, and 
theoretically, the results coincided with the reported by Benard, giving the impression that his model was 
correct. However, it is known now that Rayleigh’s theory is not adequate for explaining the convective 
mechanism observed by Benard. In fact, in Benard’s experiments, the free surface was in contact with the 
atmosphere which generates a surface tension, and Rayleigh, using a plate in the upper surface, eliminated 
the surface tension’s effects. 

It should be noted that the surface tension is not constant and it may depend on the temperature 
or contaminants in the surface. This dependence is called capillarity or Marangoni effect da ED. The 
importance of the Marangoni effect in Benard’s experiments was established by Block in [ 6 ] from an ex¬ 
perimental point of view, and by Pearson |36] from a theoretical point of view. Now is recognized that 
the Marangoni effect is the main cause of instability and convection in the Benard original experiments. 

For the foregoing reasons, we consider the physical situation of a horizontal layer of a fluid in a cubic 
container of height d (a: 3 -coordinate), of length Li (a; 2 -coordinate) and width li (a;i-coordinate). The bot¬ 
tom surface of the container and the lateral walls are rigid and the upper surface is open to the atmosphere. 
In order to describe the system, we use the Oberbeck-Boussinesq approximation |D, which assumes that 
the thermodynamical coefficients are constant, except in the case of the density in the buoyancy term, 
which is considered as being pq[ 1 — a{6 — Oa)]. Here po is the mean density, 0a is the temperature of 
the environment and a is the thermal expansion coefficient, which is positive for most liquids. More¬ 
over, we assume that the surface tension is a function of the temperature, and it is approximated by 
a = (To — l{0 — 0a)- Here, (Tq is the surface tension at temperature 9a^ and 7 is the ratio of change of 
the surface tension with the temperature (7 is positive for the more commonly used liquids). Also, the 
free surface is presumed not to be distorted, that is, the vertical component of the velocity in the free 
surface always will be zero. Then, we consider that the domain, in which the fluid is confined, is given 
by fl = (0, li) X (0, Li) X (0, d). However, the analysis developed in this paper allows us to consider a 
domain il with more general geometries, specifically, we can consider H = H x (0,d), being H a Lipschitz 
bounded domain of 


In stationary regime, the RBM system is given by the following coupling between the Navier-Stokes 
equations and heat equation: 

( po{u-V)u + Vp- pAu = Po [1 - a (0 - 0a)] g in H, 

< poCp{u-V)0 = KA0 in H, (1.1) 

[ div u = 0 in 0, 

where the unknowns are u(x) = (ui(x), M 2 (x), it 3 (x)) € d(x) € M and p(x) € M, which represent the 
velocity field, the temperature and the hydrostatic pressure of the fluid, respectively, at a point x G H. 
The constant Cp is the heat capacity per unit mass of the fluid, p its viscosity, K its thermal conductivity 
and the field g is the acceleration due to gravity. 


In order to express the system in adimensional form, we make the following changes of variables: 


Xi 


X2 


", •''2 


■) -^3 


X 3 
d ’ 


dui 


dU2 , dus , 0 — 0a / d^p 


", “2 


U 3 = —^ ^ = 


01f 


p = 


POVK 


K p 

where 0u = 0c — 0a with 0c the temperature at the bottom plate, n = —^ and v = —. Thus, removing 

PoCp Po 





On the Rayleigh-Benard-Marangoni system 


3 


the primes to simplify the notation, from dm we get 

( (u • V)u = Pr [(6 + RO) 63 — Vp + Au] in O, 
< (u • V)0 = A0 in n, 

I div u = 0 in n, 


( 1 . 2 ) 


with fl = (0, 1) X (0, L) X (0, 1), where I = 4 and L = Moreover, Pr = —, R = and 

h = ———. The number R is known as the Rayleigh number and it measures the effect of buoyancy; 


Kiy 

Pr is known as the Prandtl number and it represents the relationship between the speed of diffusion of 
momentum and the rate of diffusion of heat in the fluid, and 63 is the unit vector in the third direction, 
that is, 63 = ( 0 , 0 , 1 ). 


Let us denote by dO the boundary of O and let Fi := 917 n { 0:3 = 1} and Fq := 917 \ Fi. Then, the 
following boundary conditions are imposed: 


u = 0 on Fq, U 3 = 0 on Fi, 
^+^£-=0, ^ = l,2, onFi, 
+ Bd = 0 on Fi, 0 = 0c on { 0:3 = 0}, 
If =0 on ro\{a ;3 = 0 }, 


(1.3) 


where n = (ui, 712 , 713 ) is the normal vector pointing outward, M = B = ^ > 0, and h is the heat 

exchange coefficient of the surface with the atmosphere. 


The boundary conditions for the velocity in (IE2Di are no slip conditions on the rigid and free surface. 
The condition takes into account the Marangoni effect, which represents the mass transfer at an 

interface between two fluids due to a surface tension gradient. Conditions (101)3 say that on the lateral 
surfaces there is not heat flow (adiabatics), that on the free surface is allowed the heat flow, and that the 
bottom surface is maintained at temperature 0 c (isothermal). 

From the point of view of the existence of solution of RBM problem, recently in [33] was discussed a 
bifurcation problem in which, considering either the Rayleigh number or the Prandtl number as bifurcation 
parameters. By using the local bifurcation theory due to Crandall and Rabinowitz [5], the authors showed 
the existence of stationary solutions to the problem (fT^D-drni). which bifurcate from a basic state of heat 
conduction. For basic state we refer to the exact solution of the problem (ll.2l) - (ll.3l) . which is given by 

0 B 

Uf, = 0, 0h = 0c - 3:3 and pb = P 1 X 3 + p 2 xl. (1.4) 

1 + -D 

Previously to [35], in [a US 121125] were obtained numerical results on the existence of solutions that 
bifurcate of the basic stationary states, instability and patter formation problems, as well as a validation 
of initial and boundary conditions. However, from a theoretical point of view, no more results are available 
in the literature. The main difficulty in the treating the RBM problem (HDl-dOl), beyond the coupling 
between Navier-Stokes system and heat equation, are the crossed boundary conditions dOl) involving 
tangential derivatives of the temperature and normal derivatives of the velocity field; in fact, in order 
to define tangential derivatives at the boundary, intended in the trace sense, it is necessary regularity 
of the weak solutions; this fact involves the geometry of the domain in order to use elliptic regularity 
in Sobolev spaces IT^’P for the Laplace and Stokes equations, in polyhedral domains (see [TT] |T21 [H] 
for elliptic regularity results associated to Laplace equation and [ini 12 a 1331134] for elliptic regularity re¬ 
sults related to the Stokes equation; see also [35] for related problems associated to the Boussinesq system). 


From the point of view of optimal control theory, unlike to the Navier-Stokes stationary equations, 
results on boundary control problems in which the cost functional is subject to state equations governed 
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by RBM system are not known. Some optimal control results associated with the Boussinesq equations are 
known, see for example [Iliailll|22l[2ai29]; however, the mathematical formulation and the boundary 
control problem for Boussinesq equations differ from the RBM model in the type of boundary conditions, 
principally the condition which takes into account the Marangoni effect, as well as the dimension of 

the domain. More exactly, in [T] the authors studied an optimal control problem minimizing the turbulence 
caused by the heat convection. The states are given by 3D-Boussinesq equations with Neumann control on 
the temperature. In Enau, the authors analyzed optimal control problems for the 3D-Boussinesq equa¬ 
tions with Neumann and Dirichlet boundary controls. The results of |35] do not include control theory. In 
[mus] the authors analyzed optimal control problems associated to the 2D-Boussinesq equations. The 
controls considered may be of either the distributed or the Neumann type. In the author considers the 
approximation, by finite element methods, of the optimality system and derive optimal error estimates. 
The convergence of a simple gradient method is proved and some numerical results are given. In m 
the authors studied an optimal control problem for the Boussinesq equations, also in 2D, with Dirichlet 
control on the temperature. A gradient method for the solution of the discrete optimal control problem 
is presented and analyzed. Finally, the results of some computational experiments are presented. In the 
previous references, sufficient optimality conditions were not analyzed. Optimal control problems for the 
Navier-Stokes equations through the action of Dirichlet boundary conditions have been analyzed (see for 
instance in [HIIIIIID]). In some cases, numerical results, either by solving the optimality system or by 
optimization methods, have been obtained. 

In this paper, we analyze an optimal control problem for which the velocity and the temperature of 
the fluid are controlled by boundary data along portions of the boundary; the cost functional is given by 
a sum of functionals which measure, in the L^-norm, the difference between the velocity (respectively, the 
temperature) and a given prescribed velocity (respectively, a prescribed temperature). The cost functional 
also measures the vorticity of the flow. The fluid motion is constrained to satisfy the stationary system of 
RBM. The exact mathematical formulation will be given in Section 2. We will prove the solvability of the 
optimal control problem and, by using the Lagrange multiplier method, we state the first-order optimality 
conditions; we derive an optimality system and give a second-order sufficient optimality condition. More¬ 
over, we also study the uniqueness of the optimal solution. Beside to the solvability of the optimal control 
problem, we first prove the existence of weak solutions for RBM system with nonhomogeneous boundary 
data, as well as the uniqueness and regularity properties. It is worthwhile to remark that the proof of 
existence and regularity of weak solutions for RBM system is not a simple generalization of the similar 
ones to deal with Navier-Stokes or related models in fluid mechanics m- If fact, we are considering non 
homogeneous crossed boundary conditions involving tangential derivatives of the temperature and normal 
derivatives of the velocity field, which permit to deal with pointwise constrained boundary optimal control 
of Dirichlet and Neumann type. In [35) . the boundary conditions are homogeneous, and thus, boundary 
control problems are not considered. The non homogeneous boundary conditions are assumed in spaces of 
kind Hqq (T), r C 917, which are natural from the variational point of view; these space, which are used as 
control spaces, are closed subspaces of and satisfy the embeddings Ll^(r) ^ Hqq‘^(T) ^ L^(r). 

On the other hand, to define tangential derivatives at the boundary, intended in the trace sense, it is 
necessary to analyze the regularity of the weak solutions, in particular, it is required the TJ^-regularity 
for the temperature (cf. ([131)2 and Lemma 13.21 belowl. However, due the geometry of the domain, the 
regularity of the weak solutions, when non homogeneous boundary conditions are assumed, is a nontrivial 
subject. For that, we adapt regularity results for the Laplace equation with Dirichlet-Neumann boundary 
homogeneous conditions in corner domains of miiiKis], and some ideas of [23] to treat the Robin and 
Neumann nonhomogeneous boundary conditions. On the other hand, from the point of view of the control 
theory, as far as we known, our results are the first ones dealing with with pointwise constrained boundary 
optimal control of the RBM system, by using spaces (F) as the control spaces. We give necessary 
and sufficient optimality conditions which are a significant advance in the analysis of controlling these 
equations. In order to obtain necessary optimality conditions we use an approach which differs from the 
other ones in the case of 3D-Boussinesq equations (cf. [Il[2])- In fact, in order to derive the optimality 
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conditions, in [5] the author used a theorem of Ioffe and Tikhomorov and also he assumed a property, 
called Property C, whereas that in [1], the authors used a penalization method because in that case the 
relation control-state is multivalued. It is worthwhile to remark that in the previous references related to 
convection problems, sufficient optimality conditions were not analyzed. Finally, from the point of view 
of numerical results, since the analysis of the control problem yields variational inequalities as optimality 
conditions, the numerical analysis offers new challenges, for instance, the applicability of the semi-smooth 
Newton method in order to obtain a numerical solution (cf. [14] for numerical results in the context of 
Navier-Stokes model). 

The outline of this paper is as follows: In Section 2, we give a precise definition of the optimal control 
problem to be studied and, in Section 3, we prove the existence and uniqueness of weak solutions, as well 
as we show regularity properties. In Section 4, we prove the existence of the optimal solution. In section 
5, we obtain the first-order optimality conditions, and by using the Lagrange multipliers method we derive 
an optimality system. In Section 6 , we give a second-order sufficient optimality condition. In Section 7, 
we establish a result of uniqueness of the optimal solution. 

2 Statement of the boundary control problem 

Throughout this paper we use the Sobolev space 77™(fl), and I < p < oo, with the usual 

notations for norms || • and || ■ ||LP(n) respectively. If 77 is a Hilbert space we denote its inner 

product by in particular, the L^(r2)-inner product will be represented by If 77 is a 

general Banach space, its topological dual will be denoted by X' and the duality product by {■,-)x',x- 
Corresponding Sobolev spaces of vector valued functions will be denoted by H^(r2), H^(H), L^(f2), and 
so on. If r is a connected subset of the boundary dO, we define the trace space 

■= {t; G T^(r) : there exists g G H^dO), p|p = v, = 0}, 

which is a closed subspace of 77^ (T) (see [13], p- 397), where 77^ (T) is the restriction of the elements of 
775 (5H) to T. We also will use the following Banach spaces 

X := {u = (mi,M 2 ,U 3 ) G H^(H) : div u = 0, M 3 = 0 on Ti}, 

Xo := {u e H^(H) : div u = 0 , u = 0 on Tq, M 3 = 0 on Ti}, 
r := {S' G H\n) : S = 0 on {xs = 0}}, 

X := (u e X : u • n = 0 on To \ ( 0:3 = 0}}, 

Ho^"(r) := |v G Hj/"(r) : ^ V n = 0, V n = 0 on r \ {X 3 = 0}| . 

Moreover, if T is an arbitrary subset of dO, we use the notation (/, g)r to represent the integral /p fg dS. 
In the paper, the letter C will denote diverse positive constants which may change from line to line or 
even within a same line. 

In order to establish the boundary control problem, we consider the following stationary model related 
to (ll.2l) - (ll.3l) with nonhomogeneous boundary data: 

(u • V)u = Pr [{b + R6) 63 — Vp -I- Au] in H, 

(u • V)9 = A9 in H, 
div u = 0 in H, 
u = g on rj, u = u° on Tg, 

M 3 = 0 on Ti, 

^ + m£-^ 0 , 1 = 1 , 2 , onTi, 

^ + Bff = 0 on Ti, d = 4>2 on { 0:3 = 0 }, 

on To \ { 0:3 = 0 }, 




( 2 . 1 ) 
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where Fq = Fq U Fq with Fj n Fq = 0, the vector u° = {ui,U 2 ,u^) S HgQ^(FQ) is a Dirichlet condition 

for the velocity u on Fg; the field g = ( 91 , 92 , 93 ) G HQg^(Fj!)) is given and denotes a control for u on Fj; 

additionally, (pi £ (Fo\{a :3 = 0}) is a given function which denotes a Neumman control to temperature 

1 /2 

0 on Fq \ {a ;3 = 0}, and (/)2 G ({xs = 0}) is a Dirichlet control to temperature 9 on {X 3 = 0}. 

Suppose that Ui C Hgg (Fj), U2 C H^(To \ {^3 = 0}) and U3 C Hqq^({x3 = 0}) are nonempty sets, 
and 7 i, i = 1, 6 , are constants. Assume one of the following conditions: 

(i) 7 i > 0 for f = 1,2,..., 6, with 71 , 72,73 not simultaneously zero, and Ui, U 2 and U 3 are bounded 
closed convex sets; 


(ii) 7 i > 0 for z = 1,2,3, with 71 , 72,73 not simultaneously zero, 7 i > 0 for z = 4,5,6 and Ui, U 2 and 
U 3 are closed convex sets. 


We study the following constrained minimization problem on weak solutions to problem (EH), for fixed 
a uO G hJ/"(F2): 

Find [u, 9, g, (pi, P 2 ] G X x H^(0) x Ui XU 2 XU 3 such that the functional 


J[M,9,g„Pi,P2\ = ^||rot u||| 2 (n) + f ||u - Ud||i 2 (n) + 


2 V2(ro\{a:3=0}) 2 nr 


H2{Tp) 


( 2 . 2 ) 


is minimized subject to the constraint that [u, d] is a weak solution of ()2.1I) . Here G L2(D) 
. and 9d G L^(D) are given. 


3 Well-posedness and Regularity of Solutions for (l2.1ll 

In this section we analyze the existence, uniqueness and regularity of weak solution for system EH, 
which, as was said in Section [TJ it is not a simple generalization of the similar ones to deal with Navier- 
Stokes or related models in fluid mechanics |5U] . 

3.1 Weak Solutions for (12.111 

We introduce the bilinear and trilinear forms a:XxX^]R, c:XxXxX—:>]R, ai : H^(0.)xH^{0) —>■ 
M, bi : H^(Q) X X ^ R and Ci : X x x —h M, for the velocity and temperature: 

a(u, v)= / Vu-VvdD, c(u,v, z)= / [(u ■ V)v] • z dD, 

Jq Jq 

ai(9,W)= [ V9-VWdn, bi(9,v) = [ V9 ■ p-dO, 

Jq Jq 0x3 

ci(u,9,W)= [ [(u-V) 6 i]WdD. 

Jq 

Lemma 3.1. The following relations hold for c and Ci: 

c(u, V, z) = —c(u, z, v), c(u,v,v) = 0, Vu G Xq, Vv, zGH^(f2), (3.1) 

ci(u, 6 »,W) =-ci(u,W, 6 »), ci(u, 6 », 6 ») =0, Vu G Xq, V9,W £ H^O). (3.2) 

Proof: Considering that u G Xq, i.e. u = 0 on Fq, U 3 = 0 on Fi and div u = 0, and the normal vector 
n on Fi is n = (0, 0,1), we obtain that u • n = 0 on dO. Therefore, the proof follows as in Lemma 2.2 in 
[I5], p. 285. □ 
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Lemma 3.2. ) Assume that Q is a bounded domain with boundary dO, Lipschitz, and dO, = Fq U Fi 

with Fi C {x 3 = C}, C a constant. If 0 € then 

^vi + -^V 2 dS= [ V0-^dO, VvgXo. 

OXi 0X2 Jci 0 X 3 

Motivated by the formula of integration by parts and using Lemma 13.21 we obtain the following weak 
formulation of System m- 

Definition 3.3. A pair [u,0] G X x is said a weak solution of if 

Pra(u, v) + PrM6i(0, v) + c(u, u, v) = (/( 6 >), v), Vv G Xq, (3.3) 

ci(u, 0, W) + ai(0, + {B0, VF)p^ = (^ 1 , l^)r„\{, 3 =o} > G (3.4) 

u = g on Fj), u = u° on Fg and 0 = (j )2 on {X 3 = 0}, (3.5) 

where {f{0),v) = ( Pr {b + R0) €3 ■ v dO. 

Jn 

3.2 Existence of Weak Solutions 

In order to prove the existence of a solution to the problem (I531)-(l^ we reduce the problem to 
an auxiliary problem with homogeneous conditions for the velocity u on Fg and the temperature 0 on 
{x 3 = 0}. For that, we will use the Hopf Lemma (see Lemma 4.2 of [TH], p. 28). Notice that if 
u° G HgQ^(FQ) and g G HQg^(F))), then there exist u° G H^(9n) and g G H^(9n) such that 

u°lrg=u°, u°|an\rg= 0 , f u° • n = 0 , u° ■ n = 0 on Fg \ {^3 = 0 }, 

ilri=g, ilan\ri= 0 , g • n = 0 , g • n = 0 on fJ \ {xa = 0 }. 

Thus, the function u° + g G (9f7) and 



/ (u + g) • n = / g • n + 

IdQ Jr}. 




gn 


u'J • n = 0. 


ir} 


Therefore, by the Hopf Lemma, there exists a function which satisfies the conditions 

Ug G H^(H), div Ue = 0 in H, = u° + g on dO, 




(||U° 




l|g| 






|c(v,Ue,v)| < e||v||^i(fj) , VvgX, 


^0) 


where C = C{n, H) and e > 0 is a real number arbitrarily small. Notice that |rg= and Ug |ri= g- 
Moreover, proceeding as in Lemma l3.11 we can easily prove that the following relations hold: 

c(Ue,U,v) = -c(Ue,V,u), c(Ue, U, u) = 0, Vu,vGXo, 

Ci(Ue,6', IT) = -Ci(Ue,IT, 6»), Ci(Ue,6', 0) = 0, Vfl, IT G F. 

On the other hand, arguing as in [?], we can construct a function 0s G such that 


(3.6) 

(3.7) 


05 = (1)2 on {X 3 = 0}, ^ + B0S = 0 on Fi, ^ = ifi on Fg\{x 3 = 0}, 

on On 


ll^<5llL*(n) — 


ffi(n) - ^ 


[Ui 


T^(ro\{a:3=0}) 


Wh 


T^({x3 = 0}) 


)■ 


(3.8) 
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Here <5 is an arbitrarily small number and the constant C depends on S. 

Rewriting [u,0] G X x H^{0) in the form u = + u and 9 = 9s + 0 with u G Xq and 9 G Y new 

unknown functions, from Definition 13.31 we obtain the following nonlinear problem: Find [u,0] G Xq x Y 
such that 

Pra(u,v) + Pr M bi(0,v) + c(u, u, v) + c(u, u^.v) +c(ue,u,v) = (^f(9 + 9s),v'j 

-Pra{us,v) - Pr M bi{9s,v) - Vv G Xq, (3.9) 

ci(u, 9, IH) + Cl(u„ 0, W) + Cl(u, 9s, TH) + ai (9, IH) + = (</>i, bF)r„\{, 3 =o} 

-ci(ue, 9s, VF) - ai{9s, W) - {B9s, , VVF G H. (3.10) 

Here {f{9),V) is as in Definition 13.31 

In order to prove existence of a solution [u, 0] G Xq x F of (I3.9I) - (I3.10L we introduce the mapping 
F : Xq —>■ Xq defined by F(u) = u, u G Xq, such that [u, 6\ G Xq x F is the solution of the following 
linearized problem 

Pra(u,v) Y Pr Mbi{9,v) +c(u,u,v) +c(u, Ue,v) +c(ue,u,v) = (^f(9 + 9s),'>r^ 

-Pra(ue, v) - PrM&i( 6 » 5 , v) - c(ue,Ue, v), VvGXq, (3.11) 

ci(u, 9, FF) + ci(u„ 0, IF) + ci(u, 0^, FF) + ai(0, FF) + (p0, If)^^ = (</>i, hF)r^\^, 3 ^o} 

-ci{u„9s,W) - ai{9s,W) - {B9s,W)^^,'iW &Y. (3.12) 

In next lemma, we shall show that the operator F : Xq —>■ Xq is well-defined. 

Lemma 3.4. Let u G Xq and (j)i G Hi{TQ \ {ccq = 0}). Then there exists a unique weak solution 
[u, 9] G Xq X F of problem \3.11\) - h3.1‘A) . Moreover, the following estimates hold 

||u||^i(f2) < C ^l&l + (P + M) (^||0||_f/i(n) + ||^'5||ffi(n)) + p^||u||L4(a)||ue||L4(n)^ 

+C ^||VUe||i2(Q) + ^||Ue||^i(Q)^ , (3.13) 

ll^llii-l(n) < C + l|u||L4(Q)||05||i4(Q) + (||Ue||i4(0) + 1 + P) ||05||Hl(a)) ,(3.14) 

where C is a eonstant independent of u, u, (fi and 9. 

Proof: We consider the bilinear continuous mappings a : Xq x Xq —R and oi : F x F ^ R given by 

a(u, v) = Pra(u,v) + c(u, u, v) + c(u£,u, v), Vu, v G Xq, 

ai(0, IF) = ci(u, 9, IF) + ci(ue, 9, IF) + ai(0, IF) + (b9, If) , V0, IF G F. 

Consequently, we rewrite (13.1111 and (13.1211 as 

a(u,v) = (Zg-,v) , Vv G Xq, (3.15) 

ai(0,IF) = ^^i,IF^, VIFgF, (3.16) 

where 

(lg,v) = (^f(9 + 9s),v'^ - PrMbi(9 + 9s,v) - c(u, u^, v) - Pra(us,v) - c(ue,Ue, v), Vv G Xq, 
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(<^1, Vf) = (</.!, - ci(u, Os, VF) - ci(u„ Os, W^) - a^iOs, W^) - {BOs, VF)r3 , VVF G F. 

We can verify that the operator bilinear oi is continuous and coercive on Y and (j)i G Y'. Indeed, the 
continuity of fii and 4>i it follows from the Holder inequality and Sobolev embeddings. Moreover, the 
coercivity of di follows from (I3.2p . (13.71) and the following generalized Poincare inequality: 


h|U2(0) < C (^||Vu|U2(0) + ^ \u\^ , Vu G H\n), (3.17) 

where C = C{n,0,,Y) and E is an arbitrary portion of dO of positive measure (cf. Lemma 10.9 in |30j . 
p. 327; see also [H], p. 56). Therefore, by the Lax-Milgram theorem, there exists a unique 6 € Y which 
satisfies equation (I3.16P . Knowing 0 and inserting it in the equation (13.151) , by using the Holder inequality 
and Sobolev embeddings we can verify that the operator bilinear a is continuous on Xq and Ig G Xq. 
Moreover, from (13.11) . (13.61) and using the generalized Poincare inequality (13.171) we have that a is coercive. 
Therefore, by the Lax-Milgram theorem, there exists a unique u G Xq which satisfies equation (I3.15|) . 
Finally, setting v = u in p.l5|) . IF = d in p.l6l) and using the generalized Poincare inequality (13.171) . we 
easily obtain (13.131) and (13.141) . □ 

Now, using the Schauder Fixed Point Theorem, we will prove the existence of a fixed point of F which 
yields a solution of (I3.9I) - (I3.10I) . For that, we consider the ball 7?^ = {u G Xq : ||u||jji(qj < r} C Xq, 
where r is a positive constant such that 


r >C7 \b\ + {R + M) Ui 


'H2(ro\{x3=0}) 


+ (l|ue||L4(n) + 1 + B)\\9s\\m{Q) +||Vu£||L 2 (n) 


— I 

Pr' 


u, 




It follows from (13.131) - (13.141) that F{Br) C Br, provided <5 be small enough and Pr large enough. Moreover 
F is completely continuous. This follows from the next inequality 


||F(ui) — F(u 2 )||//i(n) < — U2||i4(Qy (3.18) 

and from the compact embedding of H^(H) ^ L'*(H), where 
C 

(l|u2||L4(n)+ ||Ue||L4(n)) 

+ C{M + 'J||Q2||L4(Q) + (||Ue||2,4(n) + l + i?)||05||^i(f2)j , 

and C is a constant independent of Ui and U 2 . 

Let us prove (13.181) . Let 9i € Y he the solution of equation (13.121) corresponding to Ui G Xq and set 
Ui = F{ui), for i = 1, 2. From (13.111) and (13.121) we obtain 

Pr a(ui - U 2 , V) +c(ui, Ui - U 2 , v) + c(ui - U 2 , U 2 , V) = -C(ui - U 2 , Uj, v) 

-c(ug,Ui - U 2 ,v) - PrM6i(0i - 02,v) + [ PrR{9i - 92)v3, Vv G Xq, (3.19) 

ai(0i -92,W) + (5(01 - 92), W) = -ci(ui,0i - 02, W) 

-Ci(ui - U 2 ,02, VF) - Ci(u„ 01 - 02, IF) - ci(ui - U 2 , Og, IF), VIF G F. (3.20) 

Setting IF = 0i — 02 in (j3.20l) and using (13.21) . (13.7L the Holder inequality, the continuous embedding 
77^(H) ^ L^(H) and the Poincare inequality, it is not difficult to obtain 

||V(0i — 02)||L2(n) < C ^||02||i/i(f2) + 11^1 “ u2||L4(n). (3-21) 

Now, using (I3.14|) and the Poincare inequality (I3.17L from (j3.21l) we obtain 
||^i-^2||ffi(n) < C'(||((>i||ppi (5||u2||L4(f2) + (||ue||i4(Q) + l + 77)||05||^i(n)) l|ui- U2||p4(Q). (3.22) 
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Setting V = Ui — U 2 in (lit.191) and using (Id.lj) . (I3.6L the Holder inequality, the continuous embedding 
H^(0) ^ and the Poincare inequality, we obtain 

Pr ||V(ui-U 2 )||i 2 (n)< (||u2||L4(n)+ \\u,\\l^q))\\ui- C{PrM + PrR)\\el-92\\H^Q)■ (3.23) 

Then, using the Poincare inequality, from (j3.23ll we get 

||ui - U 2 ||j:^i(Q) < C (||u2|U4(o) + ||ue||i4(0))||ui- U2||^4(n) + (M + P)||01- 02||ffi(n)^ ■ (3.24) 


Thus, (I3.18|) follows from (I3.24|) and (13.2211 . Therefore, the Schauder Theorem implies that F has a fixed 
point u = P(u). The field u, together with the corresponding function 6 — G Y solving the problem 
(13.121) for u = u, is a solution to the problem (I3.9I1 - (I3.10I) . We collect this result in the following theorem: 

Theorem 3.5. Let (t>i G H^To \ {xa = 0}), (j) 2 ^G H^o^{{x 3 = 0}), u° e Ho^^(rg) and g e Hj^^(rJ). 
Then there exists at least one solution [u, d] S X x FI^(fl) of problem id. ,VI) - id. 5\) provided Pr be large 
enough, and the following estimate holds 


u 


ffi(n) 


ll^ll 




(l|u°|| 


H2iTl) 






'H2(ro\{x3=0}) 


+ ll</>2 


'H 2 ({x 3 = 0 }) 


(3.25) 


where the constant C depends linearly of the parameters M, B and R. 


3.3 Uniqueness of the Weak Solutions 

The purpose of this section is to determine conditions on the boundary data and parameters which 
guarantee the uniqueness of the weak solution [u,0] G X x H^{Lt) to the problem p. 311 - (13.51) . For that, 
suppose that there exist [ui, 6(i], [u 2 , ^ 2 ] G X x weak solutions of system (13.311 - (13.511 . Then, defining 

u = Ui — U 2 and 9 — 9i — 62 , we obtain that [u, 9] G Xq x Y solves the system 


Pr a(u, v) -I- PrM6i(0, v)-b c(u, Ui, v) + c(u 2 , u, v) = / Pr P 0^3 dfl, Vv G Xq, (3.26) 

Jn 

Ci{u, 9uW) + Ci{u2, 9, W) + ai{9,W) + {B9,W)^^ =0, 'iW gY. (3.27) 

Proceeding as in Lemma 13.11 we can easily prove that if U 2 G X, u G Xq and 9 gY, then c(u 2 , u, u) = 0 
and Ci(u2,0, 0) = 0. Thus, setting v = u in (13.26p . IF = 0 in (13.2711 . and using the Holder inequality, 
Sobolev embeddings and the Poincare inequality (13.1711 , we deduce 


Pr||Vu|U2(n) <PrM||V0|U2(n)+C||Vu|h2(n)||Vui|h2(n)+C'PrP||V0|U2(n), (3.28) 

l|V0|U2(o) < C||Vu|h2(o)||V0i|h2(o). (3.29) 

Using (13.29P in p.281) . we find 

Pr||Vu|U 2 (o) < C(||Vui|h 2 (n) + (PrM -b PrP)||V0i|U2(n))||Vu|h2(n). 

Now, taking into account that [ui, 0i] and [u 2 ,02] are weak solutions to the problem (13.31) - p.51) . then from 
Theorem 13.51 we have that [ui,0i] and [u2,02] satisfy the estimate (13.25L which imply that 


Pr||Vu|U 2 (o) < C(Pr(M+P)-bl) 




l|Vu||L=(n): 


where the constant C depends almost linearly of the parameters M, B y R. Therefore, if the condition 
Pr-C(Pr(M + R) + l)\\\u°\\ 1 , , +||g|| i, ,, + ||<?ii|| i, ,, + ||^ 2 || 1 ,, , J > 0 (3.30) 

' ' ; ; II ii_f/2(rg) "®"_f/2(rj) (ro\{a:3=o}) ""ff2({2:3=0})] ^ ^ 


is satisfied, we conclude that ||Vu||= 0, and consequently u = 0, which implies that Ui = U 2 . 
Moreover, using this fact in (13.291) . we obtain that ||V0||2,2(n) = 0, and consequently 0 = 0, which implies 
that 01 =02. Thus we have proved the following theorem: 
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Theorem 3.6. Let (j)i e i ?2 (Fq \ {x^ = 0}), <p 2 S Hq^‘^{{x 3 = 0}) u° e Ho^^(r§) and g e Ho^^(rJ). If 
the condition \8.S(lt\) is satisfied, then the problem has a unique solution [u,0] S X x 

Moreover, the solution [u, 9] satisfies the estimate iS.25\) . 

Remark 3.7. Observe that the condition 


Pr - C{Pr{M + R) + 1) 


I" + l 


'H2(rl 


'Ml 


' H2 (ro\{x3=0}) 




'H 2 {{x 3 = 0 }) 


> 0 


is verified if either the functions g, (j)i and (j )2 are small, or if the coefficients M, R and B are small. 
In particular, for small values of M, R and B and boundary data = 0, g = 0, (^i = 0 and 4)2 = (^c, the 
basic solution [ut,0h,pb] given by is unique. 


3.4 Regularity 

In Subsection [321 was demonstrated the existence of a weak solution [u, 0] G X x H^{0) to the problem 
(IO)-(l?31): however, taking into account the tangential and normal derivatives of the temperature at the 
boundary, we need to prove that 9 G H^(0) (see Lemma [51^ . In this subsection we analyze the following 
regularity problem for the weak solution 9 G Given u G X, find 9 G such that 

{ -A0 = -(u-V)0 
^ + B9 = 0 

li=^i 

9 = 4)2 

where Li := { 0:3 = 1}, La := = 0} and r 2 := dO \ {Li U La} (see Figured]). 


in n, 
on Fi, 
on F 2 , 
on F 3 , 


(3.31) 


Fi 

Fz 


y 
y 
y 
y 

Figure 1 : Representation of dQ. 

In this subsection, we will use the following space 

^ 00 ^(r) = {r G L^(F) : there exists g G Hi{dO), g\^ = v, g\gQ\r = 0}- 
Theorem 3.8. Let 4>i G H^(T 2 ), 4^2 € and f G LP{0) with | < p < 2. Then, the system 





{ —A9 = / in O, 

11+50 = 0 onTi, 
||=(()i onT2, 

9 = 4)2 on Fa, 


(3.32) 


has a solution 9 G IF^’^(n). 
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Proof: We first convert the problem (13.1121) with Robin, Neumann and Dirichlet conditions, in a boundary 
problem with only Dirichlet and Neumann conditions. For that, we will adapt the ideas of [23], Section 
2. First, we consider the functions rjiyc^) and 9 defined by: 


77(0:3) = exp[R( 2 a ;3 



and 


9 = r]9. 


oroblem (13.321) 

is equivalent 

II 

<1 

1 

in 


ae 

r)n 

= 0 

on 

Ti, 

as 

dn 

= 4i 

on 

r2, 

0 = 

■ 4>2 

on 

Ts, 


(3.33) 


where / = —9ri" — 2 ^- 77 ' + 77 / and = rjipi. Taking into account that 4)2 G i?Qo^(F 3 ), we consider the 
function 4>2 & Hi (dfl) such that 4^2 = 4>2 on and ^2 = 0 on 9 f 7 \F 3 . By the lifting Theorem, we have 
that exists $2 G such that d) 2 |gQ = 4‘2- Considering 9 = 9 — $ 2 , it is not difficult to verify that 

problem (13.331) is equivalent to find 9 G VF^’P(n), such that 


-A9 = f 


0 = 0 


in n, 
on Fi, 
on r 2 , 
on r 3 . 


(3.34) 


where f = f + A$ 2 , 4’3 = — and 4>i = 4>i — In order to find the solution 9 of problem 

(13.341) . we decompose 9 as the sum 0 = 0i + 02 + 6 * 3 , where 0i, 02 and 03 solve respectively the following 
problems: 


II 

sc 

<1 

1 

in n. 

0 

II 

on Fi 

0 

II 

on F 2 

II 

0 

on F 3 


(3.35) 


{ -A02 = 0 
17 = 0 
l7 = <^i 
02 =0 


in n, 
on Fi, 
on F 2 , 
on F 3 , 


(3.36) 


{ —A 03 = 0 in n, 

^ = h onFi, 

ff = 0 onr2, 

03 = 0 on F 3 . 


(3.37) 


In order to prove the existence of 0 i, 02,03 G IF^’^(il), we require the following preliminary result whose 
proof follows from Theorem 1 in m (see also m)- 
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Theorem 3.9. If F € L^{0) with | < g < oo, then the weak solution to the problem 


{ -Aw = F 

i^=o 

w = 0 


in O, 

on Fi U r 2 , 
on Fa, 


belongs to the space 

Thus, by Theorem l3.91 if / G Lp{0.) with | < p < 2, then the system (13.351) has solution 9i G W^’P{0,). 
We remember that / = —6ri" — 2-^r]'+rjf + A^ 2 - Observe that as 0 G F[^{0) ^ L^(0) then G Lp'{0). 
Moreover since ^(cca) = exp[i 3 ( 2 a :3 — \x1)\, then 77 ' = B{2 — xz)ri and 77 " = —Br] + B^{2 — X 3 )‘^r]. Recalling 
that 0 < X 3 < 1, we deduce that 77 , 77 ' y 77 " belong to L^(0). Finally, as $2 G A‘l >2 G F^(0), and 

as by initial hypothesis / G LP(0) with | < p < 2, we conclude that / G LP{0,) with | < p < 2. Thus, 
the system (13.351) has solution 9i G VF^’P(O). 

On the other hand, observe that for finding 02, 03 G solutions of (13.361) and (I3.37|) respectively, 

we can not use directly Theorem 13.91 because these systems have nonhomogeneous boundary conditions. 
Therefore, for solving the problem (I3.36p . we first divide F 2 in four parts F 2 with 7 = 1,2, 3,4, as showed in 
Figured and then we decompose the solution 02 as the sum 02 = 02 + 0| + 02 + 02 , where 0| (7 = 1,2,3,4) 
solve: 


'-A0*=O inO, 

^ = <^1 onF^, 

^ = 0 onFiU(F2\Fy, 

02 = 0 on F 3 , 


where (j)\ is defined by (j)\ = cfi on F 2 , and ^^ = 0 on F 2 \F 2 , 7 = 1,2, 3,4. 




(3.38) 


Figure 2 : Division of r 2 . 


For solving problems (I3.38L we will adapt the ideas of [23], Section 2. In the case i = 1, we divide the 
boundary of F^, denoted by cIF^, as follows: dT\ = F^^ U F^^ U F^^ U F^'^ (see Figure[3|), and we construct 
a function ifi as a solution of the heat equation: 


'^ = At/^i 



•01 = 0 

•0i(a;i, 0,2:3) = 01(2:1,0,2:3) 


in F^ X (0, 00), 
on F 2 * X (0, 00), 
on F 2 ^ X (0, 00), 
on F^. 


1,2,4, 


(3.39) 
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Figure 3: Division of dT^- 

By standard methods (cf. [5D], Ch. 10) we can verify that there exists a solution -ijji G H^{0) for problem 
(13.391) . Moreover, considering the following function 


L 

/ (L^ — 1) 

dz - X 2 i’iixi,L,X 3 ), (xi,X 2 ,X 3 ) G r 2 X (0,oo), 


we can easily see that Ti satisfies the boundary conditions in (13.3811 (for f = 1). Moreover, taking into 
account that G we deduce that Ti G Additionally, as Ti G then — ATi G 

and consequently —ATi G LP{0) for | < p < 2. Thus, by Theorem 13.91 the solution Ti of the system 


in n, 

on Ti U r 2 , 
on Ts, 

belongs to In conclusion, considering 0\ = Ti -\- Ti, we obtain that 0\ G satisfies the 

system (13.381) for i = 1. Analogously, we can find solutions 6 ) 3 , and 63 in for problems (I3.38|) 

(for i = 2, 3,4) and (13.371) respectively. Thus, considering O 2 = 612+02+^2+^2 deduce that 02 G 
is a solution to the system (13.361) . Therefore, it was verified the existence of 0i, 02 , 03 G solutions 

of (13.351) . (I3.36P and (13.371) respectively, and the theorem is proven. □ 

Now, taking into account Theorem 13.81 we prove the following theorem which guarantees the existence 
of solution of problem (13.3111 . 

Theorem 3.10. Let (j)i G H^{T 2 ), (A 2 G i7Qg^(r3), u G X and 0 G H^{Vt) weak solution of system 
Then, the solution 9 belongs to the space 

Proof: First, observe that as u G X C H^(n) then using Sobolev embeddings we obtain that u G L®(n). 
Moreover, as 0 G V0 G L^{0) and consequently —(u • V)0 G Li{0). Thus, by Theorem 13.81 

we conclude that the problem (13.311) has solution 0 G Analogously, since 0 G V0 G 

1 T ^’2 (n), and consequently, using the Sobolev embedding (fl) ^ L^{Ll) we deduce that V0 G L^{Ll). 
Therefore, —(u • V)0 G LfiVl) and from Theorem 13.81 we conclude that the solution 0 of problem (13.31|) 
belongs to H'^(Ll). □ 


{ Afi = -ATi 
^ = 0 
fi =0 


Remark 3.11. Taking into account that the geometry of LI corresponds with a cube, we are able to obtain 
-regularity for the velocity u. For that, we can apply the results of -regularity for the Stokes problem 
in polyhedral domains (see 
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4 Existence of Optimal Solutions 

In this section we will prove the existence of an optimal solution for Problem (1^ . We define the set 
of admissible solutions of Problem (E^) as follows: 

Sad '■= {z = [u, 9, g, 01,02] S X X H^{0) X Ui X X ^3 such that 
J{z) < oo and the equations (13.311 - (13.51) hold}. 

Then, we have the following result: 

Theorem 4.1. Under the conditions of Theorem \3.5l if one of the conditions {i) or (ii) given in i2.2\) is 
satisfied, then the problem i2.3l) has at least one solution, that is, there exists at least a z = [u, 9, g, 0i, 02] € 
Sad such that 

J{i) = min J{z). 

Proof: From Theorem 13.51 we have that Sad is nonempty. Denote by (z^) = ([um,0m,gm) 01m > 02m]) C 
Sad, TO € N, a minimizing sequence for which lim Sf{zm) = min Sf{z). If one of the conditions (i) or 

m—J-OO 2(iSad 

(ii) is satisfied, then there exist constants Ci, C 2 and C 3 , independent of to, such that ||gm|| — ^i’ 

ll0imlLi/m ^ r niN ^ ^*2 and ||02mllnr*.r ^ C'a. Thus, from Theorem [33] we conclude that there 

ii 2 (^i o\fX 3 =U>j n 2 

exist constants C 4 and C 5 , independent of to, such that < C 4 and ||0m|liji(Q) < Cs- Therefore, 

since Ui, U2 and 1(3 are closed convex subsets of Hqq (PJ), H^{To \ {xs = 0 }) and Hqq'^({x3 = 0 }) 
respectively, we obtain z = [u, 0, g, 0i, 02 ] S H^(D) x H^{ 0 )x\Ji XU 2 XU 3 such that, for some subsequence 
of (zm)meN c Sad still denoted by (zm)meN, we have 

Um ^ u in H^(n) and strongly in L^(D), 2 < p < 6, (4-1) 

9m ^ 9 in H^{D,) and strongly in L\0), 2 < I < 6, 
gm ^ g in H^(rJ) and strongly in L^(rJ), 

01 m ^ 01 in \ { 2:3 = 0 }) and strongly in L^(ro \ {xs = 0}), 

02m ^ 02 in Hi({x 3 = 0}) and strongly in L^({x 3 = 0}). 

Since u^lpi = g^, Um|p 2 = u° and 9 m\{a: 3 =o} ~ it follows from the properties of the trace operators 
that u|pi = g, u|p 2 = u° and 9 |{ 3 , 3 = 0 }= 02 ) so, z satisfies the boundary conditions (13. 5|) . Moreover, since 
the third component of denoted by Um^ is equal to 0 on Pi for all to £ N, then from the continuity of 
the trace operator we obtain U 3 = 0 on Pi. Also, using (14.11) we obtain that div u™ ^ div u in L^(f7), and 
given that div Um = 0 for all m £ N, we conclude that div u = 0. Moreover, as u = g on Pj and u = u° on 
Pq, we obtain that u • n = 0 on Pq \ {^3 = 0}. Therefore, we conclude that u £ X. A standard procedure 
permits to pass the limit, as to goes to 00 , in the variational formulation (li01)-(nT31). and we obtain that z 
satisfies the weak formulation (13.31) - (13.51) . Consequently we have that z = [u, 9, g, 0i, 02 ] £ Sad, and then 

J(z) > inf J(z). 

Finally, recalling that the functional J is weakly lower semicontinuous on Sad, we have that 

J(z) = min J(z). 
zeSad 

□ 

Remark 4.2. Let [u;,,0f,] the basic solution to the problem il.2\) - il.3\} given by From Theorem \4. 1\ 

we can obtain the existence of controls [g,0i,02] £ Ui XU 2 XU 3 and a weak solution [u,0] £ X x H^(Ll) 
of the problem 51) . such that the functional i 2 . 2 \} is minimized if we consider and 9d = 9^, 

the basic state. 
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5 Necessary Optimality Conditions and an Optimality System 


In order to obtain first-order optimality conditions, we start by considering the following Banach spaces: 
G = X X X Ui XU 2 xUs and H = Xq x Y, with the usual inner products and norms. Moreover, 

if r is a connected subset of the boundary dfl, we define the trace spaces 

= {v G G^(r) : there exists g G {dO), g|p = u}, 

= |v e L^(r) : 3 g G H5(aG), g|p = v,y g • n = 0, g • n = 0 on T \ {xa = 0}, 53 = 0 on Ti 

which are closed subspaces of and H^/^(r), respectively. Also, let Ug defined by 

,n^/g onTj, 

® \u 0 onVl 

Then, taking into account that g G HgQ^(rQ) and u° G HQg^(rQ), we can easily prove that Ug G HQg^(ro). 

Also, we consider the following operators : G —>• Xg, ^"2 : G —>• Y', J'a : G —>• Hy^(ro) and : G —>■ 
1 /2 

He ({ 2:3 = 0}), defined at each point z := [u, 9, g, 4 >i,4>2] by: 

( (J'i(z),v) = Pra(u,v)-pPrM 6 i(d,v)-hc(u, u,v) - (/( 6 »),v), Vv G Xq, 

I {H 2 iz),W} = Cl (u,d, IT)+ai( 0 , IT)+ (i? 0 ,lT)r^-(</.!, VITgF, 

I Taiz) = u|p^-uO, 

[ P4(z) = 0|{,3=o} 

In order to simplify the notation, let us denote by M the space 

M = x;, X y' X Hi/2(ro) X py"({x3 = 0}), 

and define the operator 


P : G —> M, such that P(z) := [Pi(z), ^2(2), •A3(z), p4(z)]. 

Then the optimal control problem (12.2p is equivalent to: 

Find z := [u, 9, g, (j)i, (j)2] G G such that the functional 

y[u,6»,g,<))i,<))2] = ^||rot -p + TllgHj^i^pi) 

+ f 1 +?|k2fi 

2 V2(ro\{a:3=0}) 2 V2({x3=0}) 

is minimized subject to {F{z), [v, IT]) = [(Pi(z), v), (^2(2), IT), p3(z), p4(z)] = [ 0 , 0 , 0 , 0 ]. 


(5.1) 


5.1 Existence of Lagrange Multipliers 

In this subsection, we will prove the existence of Lagrange multipliers. For that, first we will establish 
a regularity condition for an optimal solution z = [u, 0, g, ^ 1 ,02] G Sad, as was established in [40], p. 50. 
We follows the ideas of M- We start by establishing the following two lemmas related to the Frechet 
differentiability of F and J. 

Lemma 5.1. The operator F is Frechet differentiable with respect to z = [u, 0, g, 0i, ^ 2 ] G G. Moreover, 
at an arbitrary point z = [u, 9, g, 0i, ^ 2 ] G G, the Freehet derivative operator of F with respect to z is the 
linear and bounded operator F^z) : G — >■ M such that at each point t = [hi, / 12 , r, g, t] G G, is defined by: 

{ (Piz(z)t, v) = Pra(hi, v) -h Pr M bi{h 2 ,v) + c(u, hi, v) -h c(hi,u, v) - Pri?(/i 2 , ^3)4,2(52), 

(P2z(z)t, IT) = ci(u, h2, IT) + ci(hi, 0 , IT) + ai(h 2 , IT) + {Bh2,W)^^ - {g, Il^)r3\(.3=o} - 
P3z(z)t = hi]p^ - Bir, 

P 4 z(z)t = /i 2 |{ 3 , 3 = 0 } -G 


(5.2) 
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for all [v, VF] e H, where Bi G £(Ho/^(ri), Hj^^(ro)) is defined hy 


Biv.= 


r on rj, 
0 on Fq. 


(5.3) 


Lemma 5.2. The functional J is Frechet differentiable with respect to z, = [u, 9 , g, (f>i, 4 > 2 ] S G. Moreover, 
at an arbitrary point z = S G, the Frechet derivative operator of J with respect to z is the 

linear and bounded operator 77z(z) : G —> K such that at each point t = [hi, h 2 , r, Q, t\ G G, is defined by: 


Jz{z)t = 7 i(roi u, rot hi)i:, 2 (n) + 72(11 - u^, hi)i 2 (n) + 73(0 - 9d, h 2 )L^n) + 74 (g, r) 1 




+ 75 (<^i, ({2:3=0}) ■ 


(5.4) 


In the next lemma, we will give a condition to assure that z G Sad satisfies the regular point condition 
(see [10], p. 50). Thereafter the existence of Lagrange multipliers is shown. 

Lemma 5.3. Let z = [u, 0, g, (/)i, ^ 2 ] G Sad be a feasible solution for the problem h5.1\) . If Pr is large 
enough and M, R are small enough such that 

fio ■■= min jpr - C (Pr(M + R) + ||u||j^i(n) + ||0||^i(a)) , ^ - CPr{R + M)| > 0, (5.5) 


where C is some positive constant, which only depends on the domain LI, then z satisfies the regular point 
condition. 


Proof: Given [a, b, c, d] G M, it is sufficient to show the existence of t = [hi, / 12 , r, q,t\ G G such that 

Pra(hi, v) + Pr M bi{h 2 ,^) + c(u, hi, v) + c(hi,u,v) - PrR{h 2 ,vfijL'^[n) = (a, v), Vv G Xq, (5.6) 

ci(u, h 2 , W) + ci(hi, 0, FF) + ai(/i 2 , VF) + {Bh 2 , FF)r, - {g, VF)r^\{, 3 ^o} = ( 6 , VF), VVF G T, (5.7) 

hilro = c + Si(r-g), (5.8) 
h2\{xs=o} = d +(t - ^ 2 )- (5.9) 

Setting [r, p, r] = [g, ^ 1 , (/() 2 ], we have that hi|r(, = c and /i 2 |{x 3 =o} = d. Then, proceeding as in the 
beginning of Subsection 13.21 we can prove that there exist [hf, h^] x H^{Ll) such that hflro - c and 
^ 2 1{ 3 , 3 - 0 } “ Therefore, rewriting the unknowns hi,/i 2 in the form hi = h| + hi, /i 2 = ^2 + ^2 with 
[hi, / 12 ] G H new unknown functions, from (lOll-dOll. we obtain the following linear system: 

Pra(hi,v) + PrMbi(h 2 ,'v) +c(u,hi,v) +c(hi,u,v) - Pr 77 (^ 2 ,'C 3 )L 2 (n) = (a,v), Vv G Xq, (5.10) 
ci{u,h2,W)+ci{hi,e,W) + ai(h2,W) + (Bh2,w) = {b,W),VW €Y, (5.11) 

\ / Fi 

where 


(a,v) = (a,v) -Pra(hi,v) - PrM 6 i(/i 2 , v) -c(u,hi,v) -c(hi,u,v) + PrR{h 2 ,V 3 )L 2 (n), 

{b, W) = {b, VF) - ci(u, hi IT) - ci(h^ d, VF) - a,ihi IT) - {Bhi IT)^.^ + (c^i, ■ 

In order to prove the existence of a solution for (I5.10I) - (I5.11L we will apply the Lax-Milgram theorem. 
For that, we consider the bilinear form ^4 : H x H —>■ M defined by 

A{[hi,h 2 ], [v, IT]) = Pra(hi, v) + Pr M 61 (^ 2 , v) + c(u, hi, v) + c(hi,u, v) - PrP(^ 2 , ■C 3 )l 2 (o) 

+ci(u,)^2,IF)+ci(hi,0;iT) + ai(ft 2 ,IF) + (pft 2 ,fF) , (5.12) 

\ / Fi 
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and J : HI —)• R defined by /[v, W] := (a, v) + (&, W). Thus, we rewrite (I5.1QI) - (I5.11I) as 

A([hi,h 2 ],[v,lT]) =/[v,VT]. (5.13) 

It is not difficult to prove that A{-, •) is continuous and / G H'. Now we prove the H x H-coercivity of A. 
For that, taking [v, IT] = [hi,h 2 ] in (I5.12L and using the Holder, Poincare and Young inequalities and 
Sobolev embeddings we get 

H([hi,h 2 ], [hi,^ 2 ]) =Pr||Vhi||^ 2 (o) + Pr Mbi(h 2 , hi) + c(hi,u,hi) - Pri?(h 2 ,^ 13 ) 1 , 2 ( 0 ) 

+ci(hi, d, / 12 ) + II V/i 2 ||| 2 (Q) + P||^ 2 ||i 2 (ri) 

> (Pr - CPr{M + R)- C||u||^qn) - l|Vhi||i 2 (n) 

+ ^1 — CPr{R + M) — IIV/i 2 ||^ 2 (Q) 

>/3o||[hi,h2]||^, (5.14) 

where /3o = Cmin |pr - C ^Pr{M + P) + ||u||//i(q) + ||0||^i(q)) , ^ - CPr{R + M)| > 0. Therefore, 

from (|5.13|) and (I5.14|) and the Lax-Milgram theorem we conclude the existence of [hi, / 12 ] G H solution 
of (I5.10D - (I5.11D . and consequently, we obtain that [hi, ^ 2 ] G X x H^{0) is solution of (I5.6I) - (I5.9I) . □ 

In the next theorem, we will prove the existence of Lagrange multipliers provided a local optimal 
solution z = [u, 0, g, 4 >i,(j)2] G Sad verifies the regular point condition (see Lemma [O]). 

Theorem 5.4. Let z = [u, 9, g, (j)i, 1 ^ 2 ] G Sad be a local optimal solution for the eontrol problem 15.1|) and 
assume 15.511 . Then, there exist Lagrange multipliers [Ai, A 2 , A 3 , A 4 ] G Hx (Hy^(ro))' x = 0}))' 

such that for all [hi,/i 2 , r, p, r] G X x H^{Q) x C(g) x C{$i) x €{^ 2 ) it holds: 

7i(rot u, rot hi)i 2 (n)+ 72 (u-Ud, hi)i 2 (Q)+ 73(0 - ^ 2 )L 2 (n) + 75 ((^i, e)^i 

+ 74 (g,i')^^(pi^+ 76 ((^ 2 ,T)^i^^^^^^j^-Pra(hi, Ai) -PrM 5 i(h 2 ,Ai) - c(u, hi, Ai) - c(hi, u, Ai) 

+PrR{h 2 , Ai 3 )l 2 (q) — ci(u, ^ 2 , A 2 ) — ci(hi, d, A 2 ) — ai(/i 2 , A 2 ) — {Bh 2 , A 2 )p^ + {g, '^ 2 )rQ\{x 3 =o} 

— (-^3, hi Ito —i3ir)^gy2(p^)y (5-15) 

w/iereC(g)xC((^i)xC((^ 2 ) = |[wi(g-g),W2(<?ii-<^i),W3(<(>2-<^2)], wi,a;2,W3 > 0, |g, (/)i, (() 2 ] G UiX W2xt^3|. 

Proof: From Lemma 15.31 z G Sad satisfies the regular point condition. Then, there exist Lagrange 
multipliers [Ai, A 2 , A 3 , A 4 ] G HI x (Hy^(Fo))' x = 0}))' such that 

i7z(z)h —(Piz(z)h, Ai)x^,Xo “ (.p 2 z(z)h, A 2 )v',V 

— (A3, p3z(z)h)^gy2(p^^y gy2(p^) — (A4,p4x(z)h)^^y2(|3,^^pj)y > 0, 

for all [hi, h 2 , r, p, r] G X x H^{0) x C(g) x C{^i) x C{^ 2 )- Thus, the proof of theorem follows from 

5.2 Optimality System 

In this subsection, we derive the equations that are satisfied by the Lagrange multipliers r] = [Ai, A 2 , A 3 , A 4 ] 
provided by Theorem l5.4l 
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Theorem 5 . 5 . (Adjoint Equations) Let z = [u, 0, g, c/ii, (/) 2 ] G Sad be an optimal solution for the control 
problem fOP and assume \5.5\) . Then, there exist functions (Lagrange multipliers) rj — [Ai, A 2 , A 3 , A 4 ] € 
M X (Hy2(ro))' X = 0}))' which satisfy, in a variational sense, the following adjoint equations 

to the control problem lEIP- 

PrAAi + (u • V)Ai — V^u • Ai — A2V0 = yirot{rot u) — 72(11—u^) in fl, 


AA 2 + (u • V)A 2 + PrMdiv 
div Ai = 0 m LI, 


\dx3J 


+ PrR A 13 = - 73(0 - Od) in LI, 


Xi = 0 on Fq 


Ai 3 = 0 onri, A 3 = 7 i(roi ux n) — Pr(VAi • n) on Fq, 


7 i(roi uxn) — Pr(VAi • n) = 0 on Fi, _BA 2 + ( VA 2 + PrM 


A 4 = — I VA 2 + PrM 


d\i 

dx^ 


VA 2 + PrM 


dXi 

dx^ 


dx^ 


' n on {x 3 = 0 }, A 2 = 0 on { 0:3 = 0 }, 


n = 0 on Fi, 


(5.16) 


n = 0 on Fg \ {x^ = 0} . 


Proof: From (I5.15P we obtain, taking \y,q,t\ = [0,0,0], that for all [hi,/i 2 ] € X x H^{Ll), 

71 (rot u,rot hi)i 2 (f 2 )+ 72 (ii-Ud,hi)i 2 (f 2 )+ 73(0 - dd, h 2 )L^(n) - Pr a(hi, Ai) - PrM bi{h 2 ,Xi) 

-c(u, hi, Ai) - c(hi,u, Ai) + PrP(/i 2 , Ai 3 )l 2 (o) - ci(u, / 12 , A 2 ) - ci(hi,0, A 2 ) - ai(/i 2 , A 2 ) 

— {P/ 12 , A 2 )r 3 “ {-^ 3 , hi |ro)(Hy 2 (pp)y — (A 4 , /i 2 \{x3=o})(^Hy'^({x3=o})y,Hy'^({x3=o}) ~ (^■1’^) 

Taking /12 = 0 in (|5.17p . we get 

7 i(rot u,rot hi)p 2 (f 2 )+ 72 (u-Ud, hi)p 2 (Q) - Pra(hi, Ai) - c(u,hi, Ai) 
c(hi,u, Ai) ci(hi,0, A 2 ) (A 3 , hi jpo)(To))^ (Po) Vhi G X, (5.18) 

and thus, using the Green formula, we obtain 

— 71 / rot(rot u) • hi c?n + 7 i / (rot u x n) • hi dS” + 72 / (u —u^) • hi dfl + Pr / AAi • hi dfl 

J Q J 0Q J Q J Q 


—Pr 


d\i 


hi dS' + / [(u • V)Ai] • hi dfl — / V^u • Ai • hi dLl — / X 2 V 9 ■ hi dLl 

JQ. Jq Jci 


Jan dn 

— (-^a, hi |r 3 )^gy 2 ^P^jy gy 2 (P^^ =0, Vhi eX. (5.19) 

Similarly, taking hi = 0 in (15.171) . we get 

13{S - 9d,h2)L‘2{n) - PrM 6 i(/i 2 , Ai) + PrR {h 2 , Xis)L^{n) - ci(u, / 12 , A 2 ) 

-ai(/i 2 , A 2 ) - (P/ 12 , A 2 )ri - (A 4 , h 2 \03=o})(MJ'^{{x 3 =o}))',MJ'^{{x 3 =o}) = ^^2 G dd^fl), (5.20) 

and thus, using the Green formula, we obtain 

'dXi 


dxa 


/i 2 dLl — PrM 


\0X3 


/i 2 dS + PrR / A 13/12 dLl 


73 / (d — 0d)h2 dLl + PrM / div 

J Q, Jq 

+ / [(u • V)A 2 ]/i 2 dLl+ f AX 2 h 2 dLl- f ^ h 2 dS - B [ A 2/12 dS 
Jn Jn Jan c'n dpi 

-(A 4 , /l 2 = 0, V/i 2 GdJ (11). (5.21) 

Observe that, if additionally we take the test functions hi € V in (15.181) and /12 G Hq{LI) in (15.201) . we get 
— 71 / rot(rot u) • hi dO + 72 / (u —u^) • hi dO + Pr / AAi • hi dO 


+ / [(u • V)Ai] • hi dO - / V^u • Ai • hi dO - / A 2 Vd • hi dO = 0, Vhi e V, 


(5.22) 
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73 /0 ~ Sd)h 2 dO + PrM f div ( dO + PrR [ Xi^h 2 dO 

Jn Jn \ 0 x 3 J Jq 

+ / [(u • V)A 2 ]/i 2 dn+ [ AA 2 h 2 dn = 0, V/i 2 e (5.23) 

Jq Jfl 

and therefore, since [hi,/i 2 ] S X X H^{0) is arbitrary, from (I5.19|) . (15.211) . (15.221) and (I5.23L we deduce 
that [Ai, A 2 , A 3 , A 4 ] satisfy, in a variational sense, the adjoint equations (15.161) . □ 

Finally, taking [hi, I 12 ] = [0,0] in (I5.15L we obtain 

74(g, + 1501, £’)/i^(ro\{a;3=0}) {{X3=0}) 

+ (A3,i3ir)^gy2(p^)y_gy2(p^)+ > 0, (5.24) 

for all [r, g, r] e C(g) x C0i) x C02). Taking r = g - g, and r = ^2 - ^2 in (15.241) . and 

recalling the definition of Bi given in (15.31) . we obtain 


(74g + ^3,g-g)(^l(pi)),p^i(pi) + ( 75(^1 +A 2,(/>1 -^l)^„i 


Thus, the optimality conditions are 


(H 2 (ro\{a;3=0}))'.ff 2 (ro\{a:3=0}) 
+ (76<^2 + A4,</>2 - > 0. 


(74g + -^3,g g)^^i ^pi^ 

( 75^1 + X2,(l)l - 


(7602 + A 4 , 02 — 02 ) 


(i?^({x3=0}))',H5({a;3=0}) 


> 

0 , 

Vg e Ui, 

(5.25) 

> 

0 , 

V 01 e U 2 , 

(5.26) 

> 

0 , 

V02 S IL 3 . 

(5.27) 


Therefore, the state equations described in m, the adjoint equations given in (15.161) and the optimality 
conditions obtained in (I5.25F(I5.27L form the optimality system of the optimal control problem (|5.1I) . 

Remark 5.6. Following the ideas of m we could suggest a semi-smooth Newton method applied to 
constrained boundary optimal control of the RBM system. However, due to the lack of sufficient regularity 
of the Lagrange multipliers for the pointwise control constraint in the optimality system, a direct application 
of the method to the infinite dimensional problem is not possible. Therefore, following the ideas of 
seems reasonable to apply the semi-smooth Newton method to a regularization of the original control 
problem, and finally to analyze the convergence of the regularized solutions to the optimal solution. 


6 Second Order Sufficient Condition 

In this section, we will analyze sufficient conditions for z = [u, 0, g, 0i, 02 ] S Sad be a local optimal 
solution. We will establish a coercitivity condition on the second derivative of the Lagrangian C in order 
to assure that an admissible point z is a local optimal solution. Here, we recall that 

£(z,rj) = J{z,) - (j^i(z), Ai)x^,Xo “ (.^ 2 ( 2 ), A 2 )y'.y - (-’^s, .7^3(z))(Hy2|.p^^^, gy 2 ^P^j 

- (A 4 , . 

We have that the Lagrange multiplier rj satisfies £[u_e]([z, ? 7 ])[hi,/ 12 ] = 0 for all [hi,/i 2 ] G X x 
that is, 

7 i(rot u,rot hi)p 2 (n)+ 72 (u-Ud, hi)p 2 (n)+ 73(0 - /i 2 )l 2 (o) -Pra(hi,Ai) - Pr M bi{h 2 , \i) 

-c(u, hi,Ai) -c(hi,u, Ai) +Pri?(h 2 ,Ai 3 )L 2 (Q) -ci(u, h 2 ,X 2 ) -ci(hi,0, A 2 ) -ai(/i 2 ,A 2 ) 

— (i?/i 2 , A 2 )p^ — (A 3 , hi |ro)(gy 2 (p^^y_gy 2 (p^) — (A 4 , ^2 = 0, (6.1) 
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for all [hi, /12] G X X In the next lemma we will establish a key estimate which is verified by the 

Lagrange multipliers [Ai,A 2] eH. 

Lemma 6 . 1 . Let z = [u, 0 , g, (/12] an admissible point for the constrained optimal control problem i 2 . 2 \) 
and assume h 5 . 5 \) . Then, the Lagrange multipliers [Ai, A2] G H satisfy 

ll'^illffi(n) + 11^211^1(0) < -^CiM[\i,9], (6.2) 

where 9 \ := -^||u||^i(f;) + -^||u —Ud||^2(o)+73ll^~^£i|lL2(o) andCi is apositive constant depending 
only on 12. 

Proof: From ( 16 .IL setting [hi,/i2] = [Ai,A2] G H, we have 

Pr||VAi|||2(o) + I|VA2|||2 (o) + -B||A2|||2(r^) = 71 (rot u, rot Ai)i;,2(o)+72(11-u^, Ai)i;,2(o) 

+73(^ - A2)L+n) - PrMbi(X2,\i) -c(u,Ai,Ai) -c(Ai,u,Ai) 

+Pri?(A2, Ai3)i2(o) — Ci(u, A2, A2) — ci(Ai, 0 , A2). ( 6 . 3 ) 

Then, by using the Holder, Poincare and Young inequalities and Sobolev embeddings, from ( 16 . 3 p we get 

(7.^2 (7 2 p 

Pr||VAi|||2(o) + I|VA2|||2 (o) < “ '^^1112(0) + “^11^-^11112(0) 

+CPr{R + M)(||VAi|||2(o) + IIVA2|||2(o)) + C'l|VAi|||2(o)||u||i/i(o) + <^73!!^ ~ ^dlli 2 (o) 
+C||VAi||| 2 (o)||^||ffi(o) + 2 ll^"^ 2 |li 2 (o), 
where C only depends on O. Thus, we can get 

(Pr - C (Pr(M + R) + ||u|Ui(o) + H^f^qo))) IIVAi 11^2(0) 

+ ^2 ~ CPr{R + M)^ I|VA2|||2(o) < C l|u||fl-i(o) + ^ 11 “ “ Ud|||2(o) + 73 11 ^ “ ^dlli2(o)^ • 

Then, since by hypothesis / 3 o = min|pr — c{^Pr{M + P) +||u||//i(o) + || 0 ||^i(o)) j ^ ~ CPr{R + M)| > 0 , 
using the Poincare inequality we conclude (IOI1 . □ 

Theorem 6 . 2 . Let z = [u, 0 , g, </>!, ^2] an admissible point for the constrained optimal control problem 
h 2 .‘A) and assume 15 ., 11 ) . If ^^^2 A 4 [u, 9 ] < 1 , where A := min{,^ ; “’ic? M.[n, 9 ], C\ are given in 
Lemma then there exists Kq > 0 such that 

,Czz[z,J 7 ][t,t] > ATolltllG, ( 6 . 4 ) 

for all t G ker(Pz;(z)). Conseguently, the point z is a local optimal solution. 

Proof: Let t = [hi, /12, r, g, r] G G. Then, the second derivative of the Lagrangian C, with respect to z 
at the point [z, rj] in all directions [t, t], is given by 

C^^[z,T]][t,t] = 7 i||rot hi||i 2 (o) +72||hi||^2(o) +73||/i2|li2(o) +74||r||^i^j,,^ +75||i?||^i^j,^^^^^^Pj^ 

+ 76 ||Tf 1 - 2 c(hi,hi,Ai) - 2ci(hi,/i2,A2). ( 6 . 5 ) 

({ai 3 = 0 }) 

Thus, by using the Holder and Young inequalities, we bound (| 631 ) as follows 

£zz[z,»?][t,t] > 7i||rot hi|||2(o) +72||hi||i2(o) +73||/i2||i2(0) +74||r||^i^^3^ +75||e||^i^j,^^^^^^Pj^ 

+ 76111 - 11 ^ 1 (^^^^ 0 }) “^(ll^illi^bn) + ll-^ 2 ||ffi(o))l|[hi,/i 2 ]ll 5 cxi 4 i(n)- 
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If t G ker(-Fz(z)), from (I5.2ll - (l5.;ill we obtain 

Pr o(hi,v) + Pr Mbi{h 2 ,v) +c(u,hi,v) +c(hi,u,v) - PrP(/i 2 , = 0, Vv G Xq, (6.7) 

ciiu,h2,W) + ci{hi,§,W) + ai{h2,W) + (P/12, bb)r, - (e, W^)ro\{. 3 =o} = «, G F, ( 6 . 8 ) 

hi|ro=Pir, (6.9) 

/i2|{3,3=o} = T- ( 6 . 10 ) 

Proceeding as in the proof of Lemma [5.31 we can prove that there exist [hf, / 12 ] G X x H^{0) such that 
hflro = 'SiI", ^2 1 { 3 , 3 - 0 } “ estimates in (18.811 remain true for 6s = h^, (pi = Q and (p 2 = r. 

Therefore, rewriting the unknowns hi, /i 2 in the form hi = h| + hi, /12 = /12 + /12 with [hi, / 12 ] G El new 
unknown functions, from (I6.7I1 - (I6.10I1 . we obtain 


A([hi,/i 2 ], [v,IT]) =/[v, VP], (6.11) 

where A is the bilinear form defined in ()5.12l) and / : El ^ R is defined by /[v, W] := (a, v) + (b, W) with 
(a,v) = -Pra{hl,v) - PrM biih^^v) -c(u,hi,v) -c(hi,u,v) + PrP(/i 2 , ^ 3 ) 1 , 2 ( 12 ), 

( 6 , VP) = -ci(u, hi IP) - ci(h^ 0, IP) - ai{hi VP) - {Bhi VP)^,^ + {g, fP)r„\{. 3 =o} • 

Moreover, arguing as in the proof of Lemma 15.31 from 1)6.111) we can get there exists C > 0, depending 
only on Pr, R, M, P, II [u, 0] O, such that 

||[hi,l^2]||H < ^ (ll[hi,/i2]llxxffi(f2) + ll^llff^(ro\{3,3=o})) ’ 
with /3o defined in (15.511 . and consequently. 


||[hi,/i2]||xxffi(n) — II[^ii ^2 ]||h + II[hi,/12]llxxffi(n) 




II^^ILp 


< c 

Thus, from (16.61) and (I6.12p we get 

Pzz[z,J7][t,t] > 7i||rot hi|||2(n)+72||hi||i2(0)+73||/i2||i2(n) + Y 

, 76 II ||2 


min{74,75,76K||r, .i||2 

^+ -III^ 1 >^ 2 ||Ix,, 1 ji(j 2 ) 


2 2 ({3,3=0}) 

“C'(||Ai||^i(q) + ||A2||^i(n))||[hi,/i2]|||-x^i{j2) 

'7. ^4 II^112 I TS II ^||2 I T6 ||^|i2 

- 2^^ 2 ^^^^^ffi(ro\{x3=o}) 2^^ IUi({x3=o}) 

+ (a “ C'(||Ai||^1(q) + ||A2 ||^1 (q))^ ||[hi,/l2]|||-,^j:^i(j2)- 

Therefore, by using estimate (EH) in Lemma 16.11 from (16.13)1 we have 

, 76 II 


£zz[z,r,][t,t] > fl|r||^i^^3) + fl 

+ (A-cf^A/f[u, 0 ] 


IH2({x3=0}) 


( 6 . 12 ) 


2 


yi 


' H2 (ro\{x3 = 0}) 


(6.13) 


ff2(ro\{2:3=0}) 
l/2\ 

1 Il[hi,/i2j 


2 ({3:3=0}) 


'xxffi(n)‘ 
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/r' ~ \ 1/2 / „ « \ 1/2 

Then, since by hypothesis C 0]\ < A, we deduce that T := A — C f ^A4[u, 6]\ >0, and 

consequently, 

^zz[z,r?][t,t] > min|T,y,y,y| ||t||^. 

Thus, we conclude the coercitivity condition (16.41) . Taking in particular t = [hi, /12, r, £», t] S ker(T’z(z)) 
with [r, p, r] G C{g) x C(^i) x C{(j) 2 ), we obtain that the point z is a local optimal solution (cf. [32]). □ 


7 Uniqueness of Optimal Solution 

In this section we will establish a result related to the uniqueness of the optimal solution of problem 
(12.211 . For that, suppose that there exist z^ = [ui ^ 2 ] & Sad (* = 1,2) optimal solutions of problem 

( 1 ^ . and let = [Ai, A 2 , Ag, A|] two Lagrange multipliers corresponding to the solutions z^ satisfying 
the relations (15.1811 . (15.2011 and (I5.25I1 - (I5.27I1 . Let us denote u = Ui — U 2 , 9 = §1 — § 2 , g — gi — g 2 , 
— ^ 1 , <^2 = <^2 “ 02 , Ai = AJ — Ag, A 2 = A 2 — A 2 , A 3 = Ag — Ag uud A 4 = A 4 — A 4 . Then, taking 
into account that Zg and Z 2 satisfy the equations (I3.3ll - (l3.5p . we deduce that z = [u, 0, g, ^g, (^ 2 ] satisfies 

Pr a(u, v) + Pr M6g(0, v) + c(ug, u, v) + c(u, U2, v) = j PrR9v3, Vv G Xq, (7.1) 

Ja 

ci(ug,0,lT)+ci(u,02,kF)+ai(0,lT) + (p0^1T\ =/<^g,VF\ , VIF G T, (7.2) 

\ / Ti \ / ro\{a;3=0} 

u = g on rj, u = 0 on Tq and 9 = ^2 on {x^ = 0}. (7.3) 

Proceeding as in the beginning of Subsection 13.21 we can prove that there exist [ue,^^] G X x H^{0) 
such that Ucjpi = g, UejgQ^pi = 0, 05||2,3^o} ~ 02, the estimates in (I3.8|l remain true for (pi = pi and 
02 = 02, and ||uc||4fi(n) < Cjlgjl^i^^^^. Therefore, rewriting [u,0] G X x H^{0) in the form u = u + 
and 9 = 9 + 9s, from (EUl-dTSI) we obtain that [u, 9\ G Xq x Y satisfies 

Pr o(u, v) + Pr M6g(6(, v) + c(ug, u, v) + c(u, U2, v) = / PrR9v3+ / PrR9sV3 

Jq Jq 

-Pra(u£, v) - Pr M6g(6(5,v) - c(u£,U2, v) - c(ug,Ue, v), Vv G Xq, (7.4) 

ciiui,9,W) + ci{u, 92,W) + ai{9,W) + {B9,WY = (Pi,w) 

\ / ro\{a:3=0} 

-ci{u„92,W) - ciiui,9s,W) - aii9s,W) - {B9s,W)^^ VVF G F. (7.5) 

Setting V = u in (17.411 . W = 9 in (17.511 . and using the Holder inequality, Sobolev embeddings and the 
Poincare inequality, we get 

Pr||Vu||i 2 (n) < C (||Vu||i2(o)||u2||4i-i(n) + {Pr + ||ui||/i-i(f 2 ) + ||u2||_f4i(n)) ||ug|| 4 i-i(o)) 

+C(PrP + PrM) (||V0|U2(a) + , (7.6) 

||V 6 (||l 2 (q) < C ^||0i||^i + (||Vu||i 2 (Q) + ||ue||g4i(n)) ||02||ii-i(n)) 

TC (1 + ||ug||4g-i(f2) + P) ||05||/4i(f2). (7.7) 

From inequality (13.2511 we have that + ||^i||i4i(n) < CSi, / = 1,2, where 


Si C^llu + llgill^i(pi) + ll0ill^i(p^\{,^3^Q}) 


( 7 . 8 ) 
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Thus, from (USD, dlZl) and (m we obtain 

Pr||Vu|U 2 (n) < CiiS 2 + PriR + M),52)||Vu|U 2 (o) + (Pr + Si + S 2 + Pr(P + M),52)||u,||^i(o)) 
+CPr{R + M)(l + S'! + P)|| 05 ||^i(n) + CPr{R + 

< C{{S 2 + Pr{R + M)^ 2 )||Vu|U 2 (n) + (Pr + 51 + ^2 + Pr{R + M)S 2 )m 


+CPriR + M)(l + 5r + 

Taking Pr large enough and M, R small enough, from dZH) we get 


||Vu|U 2 (n) < Po (||g||,,i(pj^ + 


ll-^i 


'ff5(ro\{2:3=0}) 

where Pq = C(Pr + + 52 + P + 1)(1 + Pr(P + M))/Pr, and therefore 

l|u||/i-i(n) < l|u||/i-i(n) + ||ue||H-i(n) < C'(||Vu|| 

L 2 (n) 


)> 


(7.9) 


(7.10) 


l|g| 




< 


C{no + 1 ) (llg|l,,i(ri) + MHi({.,=o}) + ll'^illffi(ro\{. 3 =o})) ' 


In the same spirit, from (EZD and (17.101) we can obtain 


I|V0||l2(o) < Pi (l|g|l,,^(p3) + + ll^illffi(ro\{:.3=0})2 

where Pi = C{Si + ^2 + P + 1 + S2'Ho), and thus 


PllffHn) < CiUi + 1) (I|g||^i(p3) + + ll^illffi(ro\{x3=o}) 


)■ 


(7.12) 


On the other hand, subtracting equations (15.181) written for u^, A2,0^, Ag, i = 1, 2, we have 

71 (rot u,rot hi)p2(o)+72(11, hi)p2(n) - Pra(hi, Ai) - c(ui,hi, Ai) - c(u, hi, Ag) - c(hi,ui, Ai) 

—c(hi, u, Ag) — Cl (hi, 01, A2) — ci(hi, 0, A2) — (Ag, hi |ro)(Hy2^p^jy Hy^(ro) “ ^ 

Taking hi = u in (17.131) we obtain 

- Pra(u, Ai) - c(ui,u, Ai) - 2c(u,u, Ag) - c(u, Ug, Ag) - ci(u,0i, A2) 

Cl (u, 0, A2) (A3, U Ipq ) (Hy^(ro))',Hy^(ro) ^11 L2(f2) 72 I|u||p2(f2). (7.14) 

Now, subtracting equations (15.2011 written for u^, Ag, A2,0i, A4, * = 1,2, we get 

130, h 2 )mn) - PrM bi{h 2 , Ag) + PrR (/12 , Ai 3 )l 2 (o) - ci(ug, /12, A2) - ci(u, /12, A2) 

-ai(/l2, A2) - (PI12, A2)ri - (A4, h2 \{X3=0}) = 0, V/12 £ P^(0). (7.15) 

Taking /i2 = 0 in (17.1511 we have 

— ag(0, A2) — (P0, A2)ri — {^‘^0\{x3=o}){{x 3 =o}))'{{x 3 =o}) ~ ci(ui,0, A2) 

-ci(u, 0, A^) - PrM 5i(0, Ag) + PrR (0, Ai3)p2(o) = -7311^11^2(0). (7.16) 

Moreover, setting v = Ag in (EID and IT = A2 in dUl), we obtain 


Pra(u, Ag)+ PrM5i(0, Ag)+ c(ug,u, Ag)+ c(u,U 2 ,Ag) = / PrP0Ai3, 


(7.17) 
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ci(ui, 0, A2) + ci(u, 02, A2) + ai(0, A2) + ( B9, A2 ) — (4>i, ^2 

/ Fi \ /ro\{a;3=0} 


(7.18) 


Adding (I7.14I1 . (I7.16L (I7.17I1 . (I7.18|) . and using (17.31) we have 

2c(u,u, A^) + c(u, u, Ai) + ci(u, 0, A2) + 2ci(u, 0, A2) 

= 7i||rot u|||2(n) +72||u||i2(n) +73||0|li2(n)- (7.19) 


Considering the optimality condition (15.251) {74& +^ * = 1>2, taking g = gi 

at i = 2, g = g2 at i = 1, and adding these relations, we obtain 

(•^3,g|ri)(^i^pijy (rj)' (7.20) 

Analogously, from the optimality conditions (I5.26I) - (I5.27I) we get 

^^’^^)(ffi(ro\{a:3=0}))',i?^(ro\{3:3=0}) “ ^ (ro\{a;3=0}) ’ ('7-21) 

{A4,^2|{.3=0})(^i({^3^3j)y < -761102115^1 (7.22) 

By using (I7.11L(|7.12I) and (16.21) we can bound the four terms of the left hand side of (17.191) as follows 


2c(u,u, A?) < ano + l)2^(X[ui,0i])i/2 ^11 

Pn ^ 




+ Il 0 i"^ 


ff^(ro\{a:3=0}) 


,(7.23) 


c(u,u,Ai) < C(Ho + l)^-^((Al[ui,0i])^/^ + (7 W[u 2,02])^7^) 

X fllglp 1 +Il02|pi +Il0l||^i 

2 ({3)3 = 0}) "^^"ff2(ro\{x3=0}) 


(7.24) 


ci(u,0,A2) < c{no + i){ni + i)Pj.{{M[vLi,ei]Y/^+ {M[n2,92]Y/‘^) 


1/2' 

0 


X 1 I|g|l^l(p3) 


ll^^ia^({.3^0}) + ll^^l 


ff^(ro\{x3=o}) 


(7.25) 


2ci(u, 0,A2) < C(i«o + l)((^i + l)^(A4[u2,02])'72 

00 

X (wkf 1 , +Il02f 1 +Il0lf 1 

l^"®"H2(rl) "^^"ff2({x3=0}) 

From (17.191) and taking into account estimates (I7.2()F(I7.26I) we obtain 

7i||rot u|||2(n) +72||u||i2(n) +73||0|li2(n) < (2:- min{74,75,76}) 

II02|0 1 + II01I 


X g 1 
' "^"ff 2 (ri) 


(7.26) 


(7.27) 


H2({3:3 = 0}) "ff2(ro\{3:3=0})^ 

where I := C'max{('Ho + 1)^, (Ho + l)(Hi + l)}-^((A4[ui, 0i])^7^ + (A4[u2, 02])^^^). By assuming Pr 

^0 __ 

large enough and R, M small enough such that Ho be small enough and I < min{74,75,70}, from (I7.27P 

we obtain 

l|u|U2(0) + ||0|U2(n) + l|g|l^i(pi) + + Il0lll44i(p„\{,3^o}) = 

which implies that Zi = Z2. Thus, we have proved the following theorem: 
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Theorem 7.1. If Pr is large enough and R,M are small enough, then the optimal solution of problem 
i2.S\} is unigue. 
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